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Dynamic second-order nonlinear susceptibilities, {2(0, co, to) = (to), are calculated here within a fully 

first-principles scheme for monolayered molybdenum dichalcogenides, 2 H-M 0 X 2 (X=S,Se,Te). The absolute 
values of across the three chalcogens critically depend on the band gap energies upon uniform strain, 

yielding the highest ~ 140 pmA^ for MoTe 2 in the static limit. Under this uniform in-plane stress, 

2 //-M 0 X 2 can undergo direct-to-indirect transition of band gaps, which in turn substantially affects x^^\o^)- 
The tunability of I’Y either compressive or tensile strainis demonstrated especially for two important 

experimental wavelengths, 1064 nm and 800 nm, where resonantly enhanced non-linear effects can be exploited: 
X^'^^ of MoSe 2 and MoTe 2 approach ~800 pmA' with -2% strain at 1064 nm. 

PACS numbers: 75.30.Gw, 75.50.Cc, 75.70.Tj 


Graphene, two-dimensional material, has attracted great at¬ 
tention for their intriguing physics such as exteremely high 
mobility ( 11 - 3 ]. However, their semi-metallic character due to 
the absence of a band gap, hinders graphene to replace for 
silicon based technology. There have been a lot of research 
efforts to artificially create band gap in graphene by introduc¬ 
ing defects, ripples, and so forth. 

M0S2, another two-dimenionsional material, has also 
revealed fascinating proprieties found in graphene: 
high mobility and mechanical strength jjsl, and even 
superconductivity | 3 . Despite many similarities to graphene, 
the hetero-atomic constitution of M0S2 naturally breaks 
valley or sublattice degeneracy. As a result, the M0S2 band 
gap is finite, as large as 1.3 eV| 113 - Other materials in 
the same family, the so-called transition-metal dichalco¬ 
genides (TMDs), or 2H-M0X2 (X=S,Se,Te), possess similar 
properties due to their equivalent structures. 

As graphene has been easily isolated from graphite, few- 
layer M0X2 has also been successfully exfoliated from bulk 
M0X2. Interestingly, while bulk M0S2 has an indirect band 
gap, the isolated monolayer-MoS2 reveals a direct band 
gapdl 0 , 3 - This indirect-to-direct transition of the band 
gap is not restricted to M0S2; other TMDs also exhibit the 
indirect-to-direct transition of the band gap in their monolay¬ 
ered structure | 3 . 

Bulk M0S2 crystallizes in the hexagonal structure with cen- 
trosymmetric space group Pb^/mmc (No. 194 ). This cen- 
trosymmetricity is broken in the monolayer, which is illus¬ 
trated in Fig. [n^-c): One unit layer has one Mo, and two S 
atoms above and below Mo, each taking A and B sites, re¬ 
spectively. The Mo-X bonds are strongly covalent, and the 
sandwich layers are weakly coupled by van der Waals interac¬ 
tions. 

Among many promising properties, the feasibility of 
achieving large second-harmonic-generation (SHG) has re¬ 
called intensive attention for an immediate practical applica¬ 
tion. In the exfoliated M0S2 film, a strong second-order non¬ 
linear optical properties is observed in the odd-layer M0S2 as 
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a consequence of broken centro-symmetricity, which vanishes 
for the even-layer ones due to the centro-svmmetricitv lfl^ . To 
date, several research groups have conducted measurements 
on the second-order non-linear optical properties on M0S2 
monolayer UMl. However, their values differ by three or¬ 
ders of magnitude from one measurement to another. The 
origin of this discrepancy is still under deb ate (Til] . Never¬ 
theless, unarguably large SHG coefficient, reaching close to 
100 pm/V, makes M0S2 very promising for non-linear optics 
applications iflol - lTil] . 

In this letter, using fully first-principles calculations of 2 H- 
M0X2 (X=S,Se,Te) monolayer, the electronic structure and 
the SHG coefficients are investigated. In-plane strain in a 
uniform manner alters band gaps, which in turn changes the 
SHG coefficients. An experimental way to realize such strain 
is proposed in Fig. Hid)- The TMD is placed above piezo¬ 
electric/dielectric materials such as PZT, whose both ends 
are pinned by metals. Applied electric field will physically 
change the shape of the substrate,which in turn changes lattice 
constant of TMD. Despite qualitative similarities of the elec¬ 
tronic structure among the 2H-M0X2, the resultant frequency- 



FIG. 1. (color online) (a) Structure of the monolayer 2H-M0X2 
(X=S,Se,Te). (b) Top- and (c) side-view of the structure of the mono- 
layer of 2 //-M 0 X 2 . Red and yellow spheres denote Mo and chalco- 
gen (S,Se,and Te) atoms, respectively, (d) A schematics to induce 
desired strain experimentally. 
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FIG. 2. (a) Band gaps (in eV) and (b) Static value of second- 

harmonic generation coefficients, = 0) (in pmA^) of 2H- 

M 0 X 2 (X=S,Se,Te) for different lattice constants. Black, red, and 
blue lines denote those of M 0 S 2 , MoSe 2 , and MoTe 2 , respectively. 
Open symbols in (a) denote magnitudes of indirect band gaps. 

dependent SHG coefficients, X*'^^(ci)), show some differences 
in their spectra. The tunability of by chalcogen atoms, as 
well as the in-plane strain, will be addressed. 

First-principles calculations are carried out using the 
highly precise full-potential linearized augmented plane wave 
(FLAPW) method llfa [ 1 ^ . The local density approximation 
(LDA)||Tj] is employed for the exchange-correlation poten¬ 


tial. Spin-orbit coupling is explicitly included in a semi- 
relativistic way. Muffin-tin radii of 2.25 (Mo) and 2.00 a.u. 
(chalcogens) are chosen. Cutoffs for wave function and po¬ 
tential representations are 14.06 and 144 Ry., respectively. 
The angular momentum expansion of charge density and po¬ 
tential inside the MT sphere are done for ^ < 8 . Summa¬ 
tion in the Brillouin zone is performed using Monkhorst-Pack 
scheme with a 45 x45 x 15 k point mesh, which includes high- 
symmetric K and M points. To reduce artificial layer-layer 
interactions, adjacent layers are separated by 20 A. We em¬ 
ployed the experimental lattice constants for the present work, 
a=3.160 A (M 0 S 2 ), 3.299 A (MoSe 2 ), and 3.522 A (MoTe 2 ), 
respectively. Internal coordinates has been fully relaxed with 
force criteria 1 x 10^^ eV/A. For simplicity, band gap correc¬ 
tion beyond LDA is not employed and the exciton effects are 
ignored. 

For calculations of the second-harmonic generation (SHG) 
coefficients, we employ the formalism by Duan et a/. lll§.[l^ . 
which is an extension of method by Sipe and Ghahramani 
and Aversa and SipeyhlThis formalism has been very suc¬ 
cessful in other works Ill^[l^l 22 l - I^ . where the imaginary 
( 2 ) 

part Xabc calculated as 
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where {a,b,c) represent the cartesian {x,y,z), and ij stands 
for the valence and the conduction band, respectively, and 
/ A j denotes the intermediate virtual state- either the virtual 

electron or hole state, = 5 (p^jiPn + PnP'ji^’ where 

Pfj = {^\Va\j) is momentum matrix element with suppressed 
wave vectors for the Bloch state indices. A 45 x 45 x 15 k 
mesh is also used for summation in Eq. ([T]) using the special 
k point scheme. The corresponding real part of XabA^’^^ is 
obtained by the Kramers-Kronig transformation of the cal¬ 
culated imaginary part. Due to the hexagonal symmetry 
in monolayered 2 H-M 0 X 2 , there is only one non-vanishing 
component in the tensor in D 3 /, symmetry: = Xxxy = 

XyS = -XmEE\- 

The calculated band gaps of 2 H-M 0 X 2 are shown in 
Fig. 12 a), where the direct (indirect) band gaps are shown in 
filled (empty) symbols. Upon pressure or lattice strain, band 
gaps decrease (increase) as the lattice expands (contracts). 
Noteworthy, there are indirect band gaps for expanded (M 0 S 2 ) 
and contracted lattices (MoSe 2 and MoTe 2 ). This feature in 
connection to the band structure will be discussed later. Be¬ 
fore we proceed, we emphasize here that ± 2 % strain is fea¬ 
sible in terms of energetics point of view. [See Supplemental 
material at [http:// .] for discussion in Sec.Il i^ 


Band gaps agree well with experiment. This is acciden¬ 
tal acknowledging LDA band gap underestimation implying 
exciton binding energy is large. Other first-principles calcu¬ 
lations utilizing many-body effect have revealed that exciton 
binding energies are 1.00 eV and 0.62 eV for MoS 2 (lltl and 
MoSe? ll2^ . respectively. That for MoTe 2 is not available at 
this point. In this work, for simplicity excitons are not fully 
taken into account. Neglect of excitons in our work will not 
alter physics qualitatively. 

The SHG coefficients for the static limit, = 0) [= 

as a function of different lattice constants are plotted 
in Fig.|2b). X*'^^(0) of M 0 S 2 increases monotonically as the 
lattice expands due to the smaller band gaps, reaching 97.04 
pm/V for a +2% change. For MoSe 2 and MoTe 2 , X*'^^(0) also 
increases as the lattice expands. However, both cases have 
dips at -1% lattice constants. This is due to smaller band gap 
of -1% lattice with respect to -2%. Moreover, indirect band 
gaps of -2% lattice are much smaller than the direct ones of 
-1% lattice by 0.2 eV. Above all, the enhancement of X^^^(O) 
in most cases is well understood: lattice expansion decreases 
the band gap, which in turn increases X*'^^(0). Interestingly, 
;K^^)(0) of M 0 S 2 reaches almost 100 pm/V, whereas those of 
MoSe 2 andMoTe 2 exceed 100 pm/V with the maximum value 
as high as 170 pm/V. 
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FIG. 3. Band structure of monolayer 2 H-M 0 X 2 (X=S,Se,and Te) for 
different lattice strains: -2, 0, and +2% with respect to the experi¬ 
mental lattice constant. Red arrow in each panel indicates whether 
band gaps either direct or indirect. 


For a more quantitative analysis, the band structures of 2H- 
M 0 X 2 are plotted in Fig. |3 [See Supplemental material at 
[http:// .] for band plots of other strains in Fig.SlI lE^ Impor¬ 
tant transitions from the valence band maximum (VBM) to 
the conduction band minimum (CBM) are labeled. A denotes 
the transition from VBM (K^i) to CBM (Xd) at the K point, 
B denotes from the spin-orbit split states {K^i) to Xd, C and 
D are transitions at M and F, respectively. Near the midpoint 
between X-F in the valence band is labelled as I). Red arrows 
connect VBM to CBM: 

As mentioned previously, indirect band gaps are shown for 
expanded (contracted) lattices of M 0 S 2 (MoSe 2 and MoTe 2 ). 
In all cases, direct band gaps occur at K. On the other hand, in¬ 
direct gaps of M 0 S 2 occur from the Fyi toKd', those inMoSe 2 
and MoTe 2 occur from Ty to K^- SOC splittings of M 0 S 2 , 
MoSe 2 , and MoTe 2 are of 151, 190, and 229 meV, respec¬ 
tively, which are ener gy d ifferences of Ky\ and Ky 2 - As shown 
in previous studies 012^, VBM and CBM consist mainly of 
Mo d states with a weak contribution of chalcogen p states. In 
D 3 /, symmetry, d states have three irreducible representations: 
xzjyz, xyfx^-y^ and z^- More specifically, at the K point, the 
VBM and CBM are predominantly derived fromxy/x^-y^ and 
z^ states, respectively. On the other hand, z^ (xz/yz) domi¬ 
nates the VBM (CBM) near F point. At T, midpoint of K- 
F,while the CBM is mainly of x^-y’^jxy state, the VBM is of a 


mixture of x^-y^fxy and xzjyz. In MoSe 2 and MoTe 2 , on the 
other hand, while dominant states in Ky\ 2 , Kc, and Fc^v are 
similar to M 0 S 2 , the valence band of the midpoint, Ty, con¬ 
sists of a mixture of xz/yz and z^. We note that indirect gaps 
in M 0 S 2 occur such that Fy shifts upward with lattice contrac¬ 
tion, while in MoSe 2 and MoTe 2 Ty moves upward with lattice 
expansion. Since the d contribution of Fc in M 0 S 2 and of Ty 
of MoSe 2 and MoTe 2 are both xzjyz, contraction of lattice in 
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FIG. 4. Second harmonic coefficient (in pm/V) of 2 //-M 0 X 2 
(X=S,Se, and Te) for different lattice strains: -2, -1, 0, -l-l, and -1-2% 
with respect to the experimental lattice constant. Blue and red dotted 
vertical lines denote half of the band gaps (//) and band gaps (£). 
Two vertical lines indicate infrared frequency range, 1064 nm (1.16 
eV) and 800 nm (1.54 eV). 

The frequency-dependent dispersion of the SHG coeffi¬ 
cients, ((£>), are plotted in Fig.[3]for different strains. The 
frequency window is chosen from 0 to 2.5 eV, which includes 
the important frequency range (IR) of 1064 nm (1.16 eV) and 
800 nm (1.54 eV). [See Supplemental material atjhttp:// .l ll^ 
to refer frequency-dependent (oi) for a wider range of fre¬ 
quencies in Fig.S2.] 

For the M 0 S 2 , in the absence of strain, the large peak 
around 2.0 eV reaches 600 pm/V just above the band edge 
around 1.0 eV. Below the band edge (£), two smaller peaks 
appear at half of the band gap {TT) and 0.2 eV below the band 
gap. From the band analysis, also in qualitative a gree ment 
with other workon the SHG coefficients of MoSo ll li M, 
these two peaks are from A/2, Bj2 and C/2, D/2, respec¬ 
tively, where the division by two indicates half of the energy 
values of the corresponding transitions. The enhancement of 
X^'^Hco) is apparent for H-1 and h-2% lattice change at a fre¬ 
quency much larger than 2.0 eV. On the other hand, large 
peaks also occur at lower frequencies for contracted lattices. 
In particular, the peak around the IR frequency range is promi¬ 
nent: x^^^ ~ 600 pm/V. 

Four distinct peaks are apparent for m < 4 eV in the case 
of MoSe 2 with zero strain. Two peaks at H and at 3 eV are 
small, and the peak just below E is pronounced, while that at 
2 eV is rather broadened. For H-1 % strain, the first peak at H is 
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sharpened, while at E it is still sharp but reduced. On the other 
hand, for +2% strain, there is a dip between peaks at H and 
E, where that at the band edge becomes broadened with third 
peak gets more pronounced. The change of spectrum between 
these two stretched strains originates from the change of the 
bands; the top valence bands remains almost unaltered, while 
the conduction band changes, though a little bit, between T- 
M and K-Y. For the contracted lattices, the direct band gap 
becomes indirect, where the CBM and the iCi very close 
in energy, hence the two peaks below E tend to merge. 

For the zero strain of MoTe 2 , two peaks are clearly no¬ 
ticeable, where the first peak is just below E, while the sec¬ 
ond peak is well above E near 2 eV. For positive strains, two 
peaks slightly shift to lower energies. More specifically, for 
+2% strain, the first peak is a little smoothened between H 
and E, whereas the second peak just below 2 eV is more en¬ 
hanced. On the other hand, for negative strains, the first peak 
is split into two - one above E, the other near El with en¬ 
hancement with respect to the zero strain case. We note that 
within IR frequency range approaches 800 pmW for 
both MoSe 2 and MoTe 2 when lattice is contracted by 2% (-2% 
strain). 

In conclusion, we predict using first-principles calcula¬ 
tions giant non-linear second-harmonic susceptibilities of the 


molybdenum dichalcogenides, 2 //-M 0 X 2 (X=S,Se,Te). We 
have found that band gaps can be altered by lattice strains; 
Expansion (contraction) of the lattice leads to a decrease 
(increase) of the band gaps. As a result, static values of 
the second-harmonic susceptibilities, X^^^(O), can be altered 
by lattice strains. Frequency-dependent SHG coefficients, 
for different strains are also investigated, where the 
peak structure strongly depends on little changes in the elec¬ 
tronic structure. Thus, we have shown that large values of 
are accessible by applying lattice strains for IR frequency 
range. In our work, the exciton spectra are ignored. If the 
formation of excitons are fully taken into account, 
should be red-shifted by amount of exciton binding energies. 
Moreover, if those excitons are robust in optical spectra such 
as absorption coefficients and dielectric function, there will be 
additional peak at the lower end of the (o). 
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